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Abstract 



We study the ground-state properties of weakly frustrated Heisenberg fer- 
rimagnetic chains with nearest and next-nearest neighbor antiferromagnetic 
exchange interactions and two types of alternating sublattice spins St > S2, 
using 1/5 spin- wave expansions, density-matrix renormalization group, and 
exact-diagonalization techniques. It is argued that the zero-point spin fluctua- 
tions completely destroy the classical commensurate-incommensurate continu- 
ous transition. Instead, the long-range ferrimagnetic state disappears through 
a discontinuous transition to a singlet state at a larger value of the frustration 
parameter. In the ferrimagnetic phase we find a disorder point marking the 



1 



onset of incommensurate real-space short-range spin-spin correlations. 
PACS: 75.50.Ee, 75.10.Jm, 75.30.Kz, 75.10.-b 
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I. INTRODUCTION 



Over the last years a great deal of interest has been concentrated on the physics of quan- 
tum Heisenberg spin chains with competing exchange couplings both for half-integer and 
integer site spins (see, e.g., Refs. [0-0 and the references therein). In general, frustration 
reduces the antiferromagnetic correlations and in some cases may produce various exotic 
quantum ground states such as the dimerized state, the spin-nematic state, or some sort 
of spin-liquid states. At the same time, little is known about the role of frustration in the 
so-called mixed spin chains constructed from two or more kinds of site spins regularly dis- 
tributed on the lattice. Recently, a number of studies concerning the physics of unfrustrated 
mixed Heisenberg chains with two kinds of spins have been published 0-0. Depending on 
the periodic array, the mixed models can form massive or massless ground states and display 
a rich variety of new features f||7|]. Meanwhile, various mixed one-dimensional quantum spin 
systems have already been synthesized in the last decade ||. 

In this report we address one-dimensional spin chains containing two different alternating 
site spins Si > S2 per unitary cell and interacting via competing antiferromagnetic neighbor 
and next-neighbor couplings, J\ > and J 2 > 0, Fig.l. The Hamiltonian of the system 
reads: 

H — J\ } J Si n S 2ri+ 5 + J2 X! (Si n Si„+i + S 2n S 2n+ i) , (1) 

re, 5 n 

where the integers n number N unitary cells and the vector 5 = ±1/2 connects the S\ 
spins with the nearest-neighbor S2 spins. The size of the elementary cell is unity. In what 
follows we frequently use the notations w = Si/S 2 , S = S 2 , J = J%j J\, J\ = 1- Theoretical 
models of quantum ferrimagnetic systems with competing interactions have already been 
discussed in the literature 0. However, these models consider complicated multiple spin- 
spin interactions which are far from the real experiment. On the other hand, little is known 
about the simple frustrated Heisenberg models containing only two-spin couplings. 
In the classical limit, the ground state of Eq. Ql|) can be described by the ansatz: 
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Si n — Si [ucos(Qn) + vsin(Qn)} , i = 1, 2, (2) 

where u J_ v are unit vectors in the spin space. The classical commensurate ferrimagnetic 
state with a pitch angle between the nearest spins 9 = Q/2 = tt is stable up to the phase 
transition point 

Jc = 2(J + 1)- (3) 

For larger values of frustration, J > Jc, the stable state is a spiral state with an ordering 
wave vector given by 

COS— = ; r. 4 

2 2J(w 2 + l) V ; 

In the limit J = oo, 9 = n/2 and the system is composed of two decoupled antiferromagnetic 
chains with site spins S\ and S2. The long-range Neel order is excluded in this limit, but one 
might expect the classical result to give a guide to possible finite-range order. In the extreme 
antiferromagnetic case, w — 1, Eq. (|4j) reduces to the well-known result for homogenous 
frustated antiferromagnetic chains, and Jc reaches its upper bound 1/4. Recent studies of 
the J = system confirm the expectation that the classical ferrimagnetic state survives the 
quantum spin fluctuations [f|]. Note that the ferrimagnetic state is characterized by both 
longitudinal ferromagnetic and antiferromagnetic long-range orders 0. Thus, a special 
feature of the discussed model, Eq. (D), is the existence of an order-disorder quantum phase 
transition. The character of the latter transition should crucially depend on the values of 
quantum site spins Si and S 2 - 

II. SPIN- WAVE ANALYSIS 

In this report we consider relatively small values of the frustration parameter J. In the 
unfrustrated case, the linear spin-wave theory (LSWT) has already been successfully applied 
in a series of recent works [Q. The second-order 1/5* series for a number of quantities were 



shown |T0| to reproduce with a high precision the density matrix renormalization group 



(DMRG) results. Thus, we can expect that at least in the weakly frustrated region the 
spin-wave theory (SWT) gives a realistic picture of the low-energy physics of the system. 
However, its application to other frustrated spin systems reveals that in many cases the linear 



theory gives better qualitative results [|Tl|]. In this respect, we observe similar behavior in 
the present system, so that the use of spin-wave expansions requires some care. 

In the ferrimagnetic region LSWT predicts the existence of two types of elementary 
excitations: 

4°' 6) = 2S ( ^±^e k ± i^zM] , (5 ) 



where a k = 1— w(l— cos k), b k = w— J(l— cos k), e k = J 1 — rjf., and r\ k = 2y/w cos(fc/2) / '(a k + 
b k ). The excitations are gapless, ~ k 2 for small k. They describe ferromagnetic 
magnons in the sector with a total spin (Si — S^iV — 1. The E^ excitations are gapful, 
= A + 0(fc 2 ), and belong to the sector (Si — S2)N+ 1 (antiferromagnetic magnons). The 
existence of gapless excitations at k — reflects the continuous symmetry of the Hamiltonian 
under a global rotation of spins. On the other hand, for w > 1 the antiferromagnetic 
excitations acquire a gap due to the violated sublattice symmetry. In the weakly frustrated 
region the structure of low-lying energy levels can also be predicted from the Lieb-Mattis 



theorem flE? ]: For every finite N, the energy levels order according to the rule: 

E(S t + l)>E(S t ), ,S t >S g , (6) 



E(S t )>E(S g ), ,S t <S g , (7) 

where St is the total spin of the state and S g = (Si — S^iV is the spin of the ground state. 
Strictly speaking, the Lieb-Mattis theorem does not work in the frustrated system, but 
due to the continuity principle we can expect that the order of states survives up to some 
finite J |n| . The above picture is also confirmed in our exact-diagonalization studies of small 
clusters. On the other hand, the strong frustration might cause a substantial disarrangement 
of levels, as is shown numerically below. 
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Exactly at the classical transition point J = Jc, the quadratic term in the small k 
expansion of vanishes, so that the spectrum reads: = ck 4 + 0(k 6 ), c > 0, J = Jc- 
For J > Jc, the classical ground state is the spiral state described above, Eq.(f|). Now LSWT 
predicts zero modes at k = and at k = ±(2ir — Q). However, the classically broken SO(3) 
symmetry in the spiral state is generally expected to be restored by quantum fluctuations 
in one-dimensional systems fl4| , so that a crucial role of boson-boson interactions should be 



expected near and beyond the classical transition point. 

It is instructive to trace back the change in the gap function A (J). The first two terms 
in the 1/5* series are: 



In Fig. 2 we show the dispersion functions, Eq. (|5]), in the case (Si, S2) = (1,1/2) for 
different frustration parameters J. Near the transition point the ferromagnetic mode is 
strongly flatted, whereas the changes in E^' are modest. In addition, Fig. 3 displays a 
smooth increase of the gap A up to the classical transition point Jc, which means that 
the antiferromagnetic excitations do not play any important role in the mechanism of the 
transition. 

In Figs. 4 and 5 second-order SWT results for the ground state energy and the sublattice 
magnetization are compared to those of the DMRG and the exact diagonalization (ED) 
numerical methods. SWT values for the energy are close to the numerical ones in a large 
region of the ferrimagnetic phase. As to the magnetization, SWT is effective only in a small 
vicinity of the unfrustrated model. Such a collapse of spin-wave series has already been 
indicated in other frustrated spin models ||15|| . We are faced with an example where the 
linear theory gives a better qualitative description. 

Let us now address the spin-spin correlations. The latter will be analysed quantita- 
tively in the framework of the DMRG method below. Here the purpose is to demonstrate 
that LSWT qualitatively captures the asymptotic behavior of the finite-range spin-spin 
correlations. As an example, we consider the spin-spin correlator K^ 2 (r) = (Si n S2 n+r + 



S x n S~2 n+r )/2. In a linear spin-wave approximation we have: 

K^r) = - 2 -^Y,-cos(kr). (9) 

The asymptotic behavior of K^(r) for long distances r depends on the analytic structure of 
the dominator in Eq. ([J). It is easy to check that in the region < J < Jp, 

J ° = 2(J+rj' (10) 

the poles in Eq. (^]) are pure imaginary. Using an appropriate standard integral, on obtains 
the following two-dimensional (2D) Ornstein-Zernike type asyptotic expression: 

r 

Ki 2 {r) ~ ^i, r^oo, J < J D . (11) 
V r 

This can be generically expected for a one-dimensional quantum problem from field- 
theoretical arguments. £ is the correlation length of the short-range transverse spin fluctu- 
ations. In the limit J = 0, LSWT predicts £ = 1/lnw, as it was also obtained by Pati et 
al. || . The LSWT function £(«/), < J < Jo, compared to the DMRG results, is shown 
in Fig. 6. A well-pronounced effect of fructration is the reduction of the correlation length. 
Qualitatively, LSWT reproduces the latter tendency. The special point Jd, known as a dis- 
order point of the first kind in classical thermodynamics, marks the onset of incommensurate 
finite-range spin-spin correlations in the chain (see, e.g., Ref. [[| and references therein). In 
the region J > Jd the poles in Eq. @ move away from the imaginary axis. Due to the 
emergence of a real part, incommensurate real-space spin-spin correlations appear. Later 
on, using also precise DMRG data, we shall discuss in more detail the physics of the disorder 
point Jo- 



in. DMRG AND EXACT-DIAGONALIZATION ANALYSIS 

To obtain quantitative results, we have used DMRG and Lanczos exact-diagonalization 
methods. DMRG is particularly adapted to the problem, for it allows treating large enough 
frustrated systems. DMRG is also free from the negative sign problems characteristic of 
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Monte Carlo methods when applied to frustrated systems. We have studied three types of 
open spin chains, (Si, S 2 ) = (3/2, 1), (1, 1/2), and (3/2, 1/2), containing up to iV = 20 unit 
cells, many times the observed correlation lengths. The DMRG results for the ground-state 
energy and sub lattice magnetization are presented in Table 1. The sublattice magnetization, 
measuring the long-range antiferromagnetic order in the ferrimagnetic state, displays only a 
slight monotonic decrease and remains finite up to the transition point Jt discussed below. 

An analysis of the correlations calculated by DMRG for large values of frustration is 
difficult due to incommensurability. As to the small frustration region, the behavior of the 
transverse spin-spin correlations can clearly be extracted for distances r up to 20 ( the limit 
is due to the very small correlation length £, leading to a very fast decay to numerical 
noise). Fits clearly reveal that the planar correlations obey the 2D Ornstein-Zernike form, 
Eq. (pi), in accord to the qualitative LSWT analysis. The fit can be done over many 
orders of magnitude. A fit to a purely exponential law is much worse. The correlation 
lengths given in Table 2 are extracted from these fits. As an example, in Fig. 6 we show 
the correlation length £ vs. J in the case (1, 1/2). We recognize a typical behavior for the 
function £(J), showing the existence of a disorder point at J ~ 0.177. Recently, it has 
been argued by one of the authors 0, that in quantum spin chains at T = the presence 
of a phase transition from a commensurate to an incommensurate phase in the classical 
limit should be generically reflected in the presence of so-called disorder points, a concept 
of classical statistical mechanics. Indeed, such disorder points have already been identified 
in a variety of cases, such as frustrated spin chains, as well as in the AKLT model 0. 
We have clearly found such quantum remnants of the classical phase transition also in the 
present model, Table 3. The point Jd is characterized by the onset of incommensurate 
real-space correlations as well as by an infinite derivative of the correlation length dC, / dJ at 
J = Jd — 0. We were not able to extract conclusively at the disorder point the crossover 
from two-dimensional to one- dimensional Ornstein-Zernike correlation functions, due to the 
extreme shortness of the correlation length. Generically, it can be expected that the disorder 
point shifts to smaller values of frustration for larger w, as it has been borne out by our 
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numerical results as well as by the spin- wave analysis, Eq. (PH). 

Not as satisfying data could be obtained for the (S 2 S z )-correlations: As we are working in 
a fixed S% otal 7^ environment, we have to subtract one-site expectation values (S z ). These 
values introduce substantial errors, given that the actual correlation is smaller by orders 
of magnitude than the subtracted value; in addition, the longitudinal correlation length is 
shorter than the transverse one. 

Finally, let us consider the structure of the low-energy levels around the classical transi- 
tion as obtained from the DMRG and ED simulations. Both methods seem to be effective 
tools revealing a complicated picture of low-energy level crossings. At the transition point 
Jt > Jc ( see Table 3) a first-order transition (level crossing) takes place from the ferrimag- 
netic state with S g = (Si — S^iV to a singlet S g = state. This is accompanied by a basic 
rearrangement of the spin correlations. Beyond the transition point Jt, the ED numerical 
analysis indicates additional level crossings. To be specific, we illustrate that for a periodic 
system with N = 6 cells and (Si, £2) = (1,1/2): The discontinuous transition from the 
ferrimagnetic ground state with total spin S g = (Si — S 2 )iV = 3 to the ground state with 
S g = appears at J = 0.229. This singlet ground state is stable up to J = 0.260. In the 
next interval (0.260,0.303) the lowest energy level is a triplet, S g = 1. Beyond the point 
J ~ 0.303 the ground state has S g = 0, but there are several crossings between singlets of 
different symmetries at J ~ 0.417, 0.735, and 0.759. These may however be due to incom- 
mensurability effects of short chains. DMRG shows that at least from J = 0.25 onwards 
the ground state has S g = 0; for values between the transition and J = 0.24, numerically 
a S g = 1 ground state cannot be clearly excluded due to a very small spacing between 
the levels just above the transition and at the same time a large DMRG error close to the 
transition. Definitely, no higher S g is observed for the ground state. Together with the ED 
result, we believe however that the ground state is indeed S g = 0. 
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IV. DISCUSSION 



In this paper the emphasis was on relatively small values of the frustration parameter 
J, when the classical ferrimagnetic state survives the zero-point spin fluctuations. As a 
whole, frustration produces a strong reduction of the finite-range spin-spin correlation length. 
Our key results are: (i) The quantum spin fluctuations destroy the classical continuous 
transition and generate a discontinuous transition from the long-range ordered ferrimagnetic 
S g = (Si — S^iV state to a singlet state, S g = 0. The transition point Jt lies beyond the 
classical transition point Jc- (ii) The discontinuous transition point Jt is preceded by 
a disorder point Jd < Jt, marking the onset of incommensurate real-space short-range 
correlations. For smaller values of frustration, J < Jd, the transverse short-range spin-spin 
correlations show typical 2D Ornstein-Zernike behavior. 

Even in the small J region, a number of open issues can be indicated. One of the 
important questions concerns the nature of the S g = ground states established beyond 
the ferrimagnetic phase. As mentioned above, the classically broken SO (3) symmetry in 
the spiral state is generally expected to be restored by quantum fluctuations in the one- 
dimensional systems. Thus, we are enforced to look for possible magnetically disordered 



states. A valuable information can be obtained from the Lieb-Schultz-Mattis theorem 16 



adapted to mixed chain spin models ||. In our case, the theorem is applicable to the systems 
(Si, S2) = (1, 1/2) and (3/2, 1), and says that the model (|l|) either has gapless excitations 
or else has degenerate ground states. Thus, we should look for phases with presumably 
some discrete symmetry broken. An analysis of possible quantum spin phases has already 
been implemented in the context of the model of frustrated ferromagnetic chains (Ji < 0, 
J2 > 0) fl7| , |l8fl . In particular, two possible nematiclike phases with additional symmetry 
breaking of reflections about a bond or about a site have been suggested |L7j . It was argued 
above that the antiferromagnetic excitations in our model do not play any important role 
near Jt due to the large energy gap A near the transition. Therefore, the instability of the 
ferrimagnetic state is expected to capture the main features of the ferromagnetic instability 
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appearing in the frustrated ferromagnetic chain. However, in contrast to the cascade of 
transitions predicted for the ferromagnetic chain ||18|| , our numerical results point towards 
a direct transition from the ferrimagnetic S g = (S± — S^iV state to a state with S g = 0. 
Such a difference in the behavior of both models at Jt may be attributed to the low-energy 
multimagnon bound states appearing in the ferromagnetic chain. In addition, the variety of 
possible nonmagnetic quantum phases beyond Jt may be considerably enlarged due to the 
existence of two kinds of site spins, Si and S2, in the mixed chain model. The complicated 
picture of low-energy level crossings described above points towards a rich phase diagram of 
the discussed ferrimagnetic model, as well. Studies in this direction are in progress at the 
moment. 
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FIGURES 



FIG. 1. The mixed-chain Heisenberg spin model studied in the paper. J\, J2 > 0. 

FIG. 2. The dispersion functions E^' b ^ in a first-order spin-wave approximation for different 
values of the frustration parameter J. 

FIG. 3. The antiferromagnetic gap A vs. J obtained in a first-order spin-wave approximation 
(solid curves) and by the Lanczos exact-diagonalization method. 

FIG. 4. The ground-state energy per cell E of the frustrated (1,1/2) ferrimagnetic chain 
obtained by the second-order SWT, DMRG, and the exact-diagonalization methods. 

FIG. 5. On-site magnetization of the first sublattice mi in the frustrated (1,1/2) ferrimag- 
netic chain obtained by LSWT, DMRG, and the exact-diagonalization methods. The second-order 
spin-wave results are not conclusive, as explained in the text. 

FIG. 6. The transverse spin-spin correlation length in the frustrated (1,1/2) system vs. J 
obtained from the LSWT and the DMRG method. Jd ~ 0.177 is the disorder point marking the 
onset of incommensurate real-space finite-range spin-spin correlations. 
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TABLES 

TABLE I. DMRG results for the ground-state energy per cell Eq, and the on-site sublattice 
magnetization m\ of the first sublattice in the systems (Si, 5*2) = (3/2, 1), (1, 1/2), and (3/2, 1/2). 



(1,1) (M) (§,5) 



J 


Eq 


mi 


Eq 


mi 


Eq 


mi 


0.00 


-1.93096 


1.14428 


-1.45409 


0.79249 


-1.96723 


1.35743 


0.02 


-1.90345 


1.13828 


-1.43337 


0.78995 


-1.92173 


1.35336 


0.04 


-1.87614 


1.13198 


-1.41278 


0.78730 


-1.87639 


1.35319 


0.06 


-1.84906 


1.12535 


-1.39234 


0.78455 


-1.83121 


1.35092 


0.08 


-1.82221 


1.11838 


-1.37206 


0.78167 


-1.78621 


1.34853 


0.10 


-1.79561 


1.11099 


-1.35194 


0.77866 


-1.74140 


1.34604 


0.12 


-1.76928 


1.10319 


-1.33199 


0.77551 


-1.69679 


1.34339 


0.14 


-1.74325 


1.09489 


-1.31223 


0.77272 


-1.65239 


1.34061 


0.16 


-1.71755 


1.08608 


-1.29267 


0.76871 


-1.60822 


1.33767 


0.18 


-1.69219 


1.07668 


-1.27332 


0.76504 


-1.58472 




0.20 


-1.66722 


1.06662 


-1.25420 


0.76115 






0.22 


-1.64267 


1.05582 


-1.23533 


0.75702 






0.24 


-1.61859 


1.04418 


-1.220 








0.26 


-1.59502 


1.03157 


-1.214 








0.28 


-1.5720 


1.01785 


-1.212 









15 



TABLE II. DMRG results for the transverse spin-spin correlation length £ in the systems 
(Si, S2) = (3/2, 1), (1, 1/2), and (3/2, 1/2) for different values of the frustration parameter J. 



J 


0.00 


0.02 


0.04 


0.06 


0.08 


0.10 


0.12 


0.14 


0.16 


0.18 


0.20 


(3/2,1) 


1.74 


1.67 


1.60 


1.53 


1.45 


1.37 


1.29 


1.20 


1.11 


1.01 


0.89 


(1,1/2) 


1.01 


0.97 


0.93 


0.89 


0.84 


0.80 


0.74 


0.69 


0.61 


0.46 


0.50 


;3/2,l/2) 


0.75 


0.72 


0.67 


0.64 


0.56 


0.53 


0.45 











TABLE III. DMRG results for the disorder point Jd and the first-order transition point J? 
found by DMRG in the systems (Si, S 2 ) = (3/2, 1), (1, 1/2), and (3/2, 1/2). J c = w/2{w 2 + 1) is 
classical commensurate-incommensurate transition point. Jg" = l/2{w + 1) is the LSWT result 
for the disorder point. 

(Si,S 2 ) J£T J D {DMRG) J c Jt(DMRG) 

(3/2,1) 1/5 0.235(5) 3/13 0.280(1) 

(1,1/2) 1/6 0.177(2) 1/5 0.231(1) 

(3/2,1/2) 1/8 0.14(2) 3/20 0.161(1) 
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Fig. 1: N. B. Ivanov et al. 



Fig. 2: N. B. Ivanov et al. 




Fig. 3: N. B. Ivanov et al. 
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Fig. 4: N. B. Ivanov et al. 
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Fig. 6: N. B. Ivanov et al. 
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